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Formal derivation of criterion for selection of superfluid phases of He in 
aerogel is presented. At the strength of the derived criterion variation of the 
order parameter of B -like phase in magnetic field differs from that of B-phase 
of pure 3 He. Possible observable consequence of this difference is discussed. 
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1. INTRODUCTION 

In pure 3 He two superfluid phases - A and B are observed. Existence 
of more then one superfluid phase is a manifestation of unconventional char- 
acter of Cooper pairing in this liquid. In both phases Cooper pairs have 
angular momentum I = 1 and spin s = 1. Their order parameters are 3x3 
complex matrices A„j, index \x refers to spin and j to the orbital projection. 
The particular form of order parameter for the A-phase is 

A£f M = A-±=d fl (m j +ih j ), (1) 

where m, n are mutually orthogonal unit vectors in orbital space and d - a 
unit vector in spin space, A is an overall amplitude. For B-phase the order 
parameter is 

A™ = A-Le i *R lxj , (2) 

where R„j is a real orthogonal matrix. In a presence of aerogel the prop- 
erty of superfluidity survives ^ and two superfluid phases (called A-like and 
B-like) are observed as well. In a simplified form aerogel can be repre- 
sented as an array of randomly distributed and randomly oriented strands 
of diameter d ps 30A. For 98% aerogel the average distance between the 



strands is I ~ 200A, i.e. comparable with the correlation length of the con- 
densate which varies depending on a pressure in the interval 160-500A. 
Aerogel makes condensate of Cooper pairs spatially nonuniform. Scattering 
of quasiparticles on the strands has depairing effect on components of the 
order parameter which extends over a distance of the order of £o ^- An array 
of strands induces random variation of the order parameter in space. It is 
convenient to represent the order parameter Aaj (r) as a sum of its average 
value and a fluctuating increment: A^j{r) = A^j + a w -(r). By the definition 
the average o^j(r) =< o^j(r) >=0. Averaging is assumed over a volume 
with a characteristic size exceeding a length for which the distribution of the 
strands is correlated. Onset of the long-range order is formally expressed 
by the condition A^j j^O. The form of A^j describes the change of symme- 
try at the phase transition and symmetry properties of the emerging phase. 
Identification of A^j for A-like and B-like phases is a key question for theo- 
retical interpretation of their properties. Fluctuating part o«j(r) influences 
averaged properties of the superfluid phases as well and because of nonlinear 
dependence of free energy of these phases on the order parameter fluctuations 
can influence a choice of A^j . This influence has been discussed recently^EI. 
In the present paper some of the previous discussion is summarised and new 
application of the developed approach is presented. 



2. RANDOM FIELD 

Since the average distance between the strands is of the order of £o per- 
turbation of the order parameter at a given point appears as effect of several 
neighbouring strands. In a vicinity of T c correlation length £(T) S> £o and 
effect of aerogel can be represented by a local random field. Its form is de- 
termined by symmetry argument and natural assumptions about properties 
of aerogel. Most important assumption is that the scattered quasiparticles 
do not exchange their spins with the strands. In practice this condition is 
secured by a coverage of the strands by one-two monolayers of atoms of 4 He. 
With this precaution aerogel can influence directly only orbital part of the 
order parameter. A principal contribution of interaction with aerogel to the 
free energy functional is of the second order in the order parameter and it 
can be written in a form 

F v = N(0) j nj^A^A^r, (3) 

where iV(0) is the density of states on Fermi level and r]ji(r) - a random 
static tensor field. Further properties of r]ji(r) can be found on a strength of 
general argument. Since F v must be real rjji(r) is Hermitian matrix and to 



secure t — > — t invariance of the energy it must be symmetrical. So ^(r) is 
real and symmetrical tensor field. Characteristic amplitude of the random 
field can be estimated with the use of results of the paper 4 . By the order 
of magnitude it is \r)ji\ ~ {^d/l 2 ). For 98% aerogel this ratio is ~ 1/10 
and effect of aerogel can be treated as a perturbation. With the account of 
interaction term (3) Ginzburg-Landau (GL) functional takes the following 
form: 

5 

Fgl = N(fl) [ d^rA^A^ + rtji^A^ + ~ £ (3J S + 

J s=l 

2 \ dxj dxj dxj dxi dxj dxi J J ' 

where r = (T — T c ) /T c , I s - 4-th order invariants in the expansion of the free 
energy over A^j (cfP). Coefficients Pi, ...0$, Ki, K2, K% are phenomenologi- 
cal constants. In what follows for their evaluation the weak coupling values 
are assumed. Possible random corrections to the gradient terms contain ex- 
tra power of the wave-vector and will not influence the following discussion. 
Comparison of two first terms in the functional indicates a physical meaning 
of the tensor rjji(r) - its isotropic part \f] mrn {r)8ji describes local variations 
of T c = T c (r) caused by fluctuations of the density of scatterers. Anisotropic 
part i]ji(r) — \r} mm {r)5ji = r)ff describes the local splitting of T c for differ- 
ent projections of angular momenta caused by the local breaking of spherical 
symmetry by the aerogel strands. 



3. CRITERION FOR SELECTION OF THE PHASES 

The order parameters of stable phases are found as minima of the func- 
tional (4). Because of a presence of the field r]ji(r) the gradient terms can 
not be omitted. According to the procedure developed for a scalar order pa- 
rameter in a scalar random field™ the order parameter A w (r) = Anj + a w - (r) 
is substituted in the GL equation, corresponding to the functional (4) and 
terms up to the second order in a^j and rfji are kept. After averaging over 
coordinate r one arrives at the equation 



x 1 „ r 9L 1 



81 

3=1 ~~ m 

d 3 i s 

'dA^dAtndAfl 



d 3 L 



dA*,, 2^dA* Kj dA vn dApi 



< 0*vn a l3l >)]=-< a ^Vlj > 



(5) 



Along with A^j the obtained equation contains averages of binary products 
of fluctuations < a un api > etc.. Equations for a^j are obtained by separation 
of fast varying terms in the GL equation 



T(h " + oX^4 ~ — — ""» + - ■' — " 



"■' r 2 ^ ys l dA*,dA vn " un ^ 8A*Ja* 



2 \ dxf dxidxj 

To proceed further one needs to introduce a nontrivial modification 
of the procedure of RefP caused by a more complicated form of the order 
parameter and its higher degeneracy. Eq.(6) and its complex conjugate form 
a linear inhomogeneous system. The corresponding homogeneous system has 
nontrivial solutions as a consequence of the degeneracy of the unperturbed 
GL functional (functional (4) without the random field) with respect to 
orbital rotations of the order parameter. These solutions are the increments 
of A^j and .A* ■ at the infinitesimal rotation 6 n : 

= $n e "' A^rjLU^j = Q n e? A^ r , (7) 

where e jnr is the absolutely antisymmetric tensor. Eq. (6) is solved by 
Fourier transformation. For further argument only a character of singularity 
of ciuj-(k) at k — > is important. To find it one can neglect anisotropy 

of the gradient terms in Eq. (6) and to substitute \K ( ^ ) instead 
of \k(^^- + 2S^-\. Then Eq.(6) is multiplied by w*,. The result 



is summed with its complex conjugate and the obtained equation is solved 
with respect to the projection a w (k) = a w -(k)a»*j + a* • (k)a; w - : 

^(k) = -^(^-Q r/ r / ;; ) (k)), (8) 

where Q r ; = A^A*; + A^A*^ . This projection gives a contribution to the 
averages < a un api > which is proportional to 

J < a w (k)a w (-k) > d 3 k ~ 

e^QHBge^'Qsn J < ^(k)^(-k) > ^. (9) 

Form of the correlation function in the integrand after averaging over direc- 
tion of k is determined by the symmetry properties of rju : 



< ^(kJi/W (_ k ) > ^ = $(fc)[2%<5 mn - 3(5 lm 5 jn + <W<5 im )]. (10) 



4tt 



Finally 

J < a w (k)a w (-k) > d 3 k ~ 



0p9q[e P e 9 QrlQsn ~\~ (<5p<j^rs <5ps<5<jr )QrnQsn] / — ^2 (H) 

The integral in the right-hand side diverges at the lower limit and the aver- 
aged products < a^a^ > give diverging contributions to Eq. (5). To make 
the iteration scheme Eqs. (5), (6) consistent one has to require the coefficient 
in front of the diverging integral to be zero. Equating to zero the expression 
in square brackets in Eq. (11) one arrives at Q rn = q5 rn , where q is a real 
number. In terms of the order parameter this condition reads as: 

A ur A U j -\- A v jA ur = const • 5 r j. (1^) 

This is a necessary condition for existence of a superfluid phase with the order 
parameter A ur . The phases meeting this criterion are "quasi-isotropic" i.e. 
the energy of their interaction with aerogel does not change at the arbitrary 
orbital rotation: 

JlfAtfA^ = 0. (13) 

Tensor of superfluid densities for these phases is isotropic and as a conse- 
quence their order parameter is not oriented by a superflow. It has been 
shown before^ that when condition (13) is met the general argument of Imry 
and MeP does not apply and the long-range order is not destroyed by the 
random field rjji(r). It might be well to point out here the difference between 
the argument of Imry and Ma and the argument presented in this section. In 
the former case the order parameter is assumed to be fixed by minimization 
of free energy without the random field. Possible disordering is caused by a 
random walk of the order parameter over its space of degeneracy. No other 
change of the order parameter is admitted. In the latter case all variations 
of the order parameter Auj are considered. That makes possible an adjust- 
ment of the order parameter to the random field. To emphasize a property 
of phases, which satisfy criterion (12) or (13) to resist disordering action of 
the random field a term "robust" was suggested^. 

Now let us discuss a region of applicability of the above argument. 
Dimensionless distance from T c : \r\ = (T c — T)/T c is limited from above by 
the condition of applicability of the GL-expansion \r\ <C 1. There exist a 
limitation from below. Even when criterion (12) is met average fluctuation 
of the order parameter diverge when |r| — > as l/^/[r|. As in conventional 
superconductors^ fluctuations are small when |r| 3> (A|? or . r //f r ^o) 2 - Here 



A corr is the correlation length for the random field rjji(r), lf r a transport mean 
free path for the Fermi excitations. Assuming A 

corr 

500A, l tr ~2000A, 

£o ~200A, one arrives at inequality |r| 3>l/30. This estimation is not very 
accurate since the poorly known value of A corr enters it in the sixth power. 
More reliable estimation comes from the experimentally observed smearing 
of the specific heat jumpD-D ft renders |r| ~l/25 for the lower limit. 

The order parameter of ABM-phase Eq. (1) does not satisfy the crite- 
rion (12). According to the above argument it can not be the order parameter 
of the A-like phase. Possible forms of the order parameters suitable for a 
description of that phase were discussed elsewher J^. For the BW-order pa- 
rameter criterion (12) is met. Identification of the B-like phase as BW-phase 
agrees with the recent pulse-NMR experiments and with the observation of 
homogeneously precessing domain (HPD) in that phased Restriction im- 
posed by criterion (12) becomes essential when B-phase is perturbed by 
magnetic field. Joint effect of aerogel and magnetic field on the BW-order 
parameter is discussed in the next section. 



4. B PHASE IN MAGNETIC FIELD 

In magnetic field two leading terms in the strength of the field have to be 
added to the free energy Eq. (4): linear in magnetic field and proportional 
to imaginary part of the tensor A^A^y 

f ( H ] = Ke^xArfAtjHx, (14) 



and quadratic in the field and proportional to real part of the same tensor 
fjp = 7; K (A[ijAlj + A V jA* ,)H^H V . (15) 



Derivation of criterion (12) is based on the invariance of GL free energy 
with respect to orbital rotations of the order parameter. The added terms 
preserve this property and as a consequence the criterion (12). This criterion 
imposes restriction on possible changes of the order parameter of B-phase. 
As in the case of pure 3 He ^ let us search the order parameter in magnetic 
field in a form of B2-phase (or oblate state): 



A 



(d/i + ie,j)(mj - irij)+ 

A u(^ - ie^irrij + inj) + 2A n f ll l j (16) 



Amplitudes Ajf, A+i , An are found by minimization of the free energy. The 
order parameter (16) satisfies criterion (12) if 



Aft = \ (Aft + Aft) . (17) 

Separate analysis^ shows that the ansatz (16) does not impose additional 
restrictions on possible forms of the order parameter, provided a deviation 
from the bulk B-phase is small. Let us introduce following parametrization: 
A n + A u = 2A, A n - A u = 25, then from eq. (17) Aft = A 2 + 5 2 . 
Minimization of the free energy in a principal order on a small deviation 5 
renders: 

r ~( H (10) 

2A[3fr + ft + 1 ' 

Coefficient C, is proportional to the derivative of the density of states over 
energy, this makes 5 really small. Keeping only principal order correction 
we arrive for the order parameter to the expression: 

A^j = AR^j + i5(e^d v - d^e u )R uj . (19) 

Eq. (19) shows that in a presence of aerogel magnetic field produces a change 
of the BW-order parameter which is different from that in bulk 3 He. In the 
bulk the order parameter remains real and principal effect is suppression of 
At i pairing amplitude while A-t-t ~ ^a- in aerogel the increment of the 
order parameter is imaginary and principal change is appearance of a dif- 
ference 5 between A^j- and An. In both cases the increments are small and 
hardly observable in a moderate field, but there is an indirect manifesta- 
tion of the induced change. Rotation matrix Raj entering Eq. (2) is often 
parametrized by an angle of rotation 6 and direction of the rotation axis n. 
The angle is found by minimization of the dipole energy^ : 



Id = t9d{A^A) 3 + A M A^). (20) 



while direction of n remains arbitrary. This degeneracy is lifted when B2- 
order parameter (16) is substituted in the dipole energy (20). For the order 
parameter (19) the relevant part of the dipole energy is 

/n = (9/2)g D 5 2 n 2 z . (21) 

It has a minimum when n z = 0, i.e. when n _L H in a contrast to the pure 3 He 
where the minimum configuration is ±n || H. A strength of the orientational 
effect is also different. Convenient characteristics of this strength is a healing 
length i-e. a distance over which orientation of n recovers to its bulk 



value being perturbed by for example a boundary of a sample. The healing 
length is found by comparison of the orientational energy with the gradient 
energy. For the orientational energy in aerogel eq. (21) one arrives at: ~ 
£o /r\ ~Fu - Dimensionless dipole constant (cf. 8 ) ~ 5 ■ 1CT 7 . Coefficient 
£ can be evaluated from the splitting of Ai - A2 transition in pure 3 He as 
10~ 3 l/koe. For H~l koe and not too close to T c healing length is of 
a few millimeters. This distance is comparable with dimensions of samples 
used in most experiments. It means that orientation of n is determined by 
geometry of samples and not by the field. 



5. CONCLUSIONS 

Properties of superfluid phases of 3 He in aerogel to a great extent are 
determined by the spatial fluctuations of the order parameter induced by 
aerogel. This contrasts with the pure superfluid 3 He where only average value 
of the order parameter is important. In a presence of aerogel contribution of 
fluctuations to physical quantities remains small only for a special class of 
order parameters. Formal criterion for selection of such phases in a vicinity 
of T c follows from very general assumptions about properties of aerogel. The 
derived criterion rules out ABM-order parameter as a possible candidate for 
A-like phase and imposes restriction on a possible variation of BW-order 
parameter in magnetic field. Because of generality of the proposed scheme 
it may be applied for description of the influence of a quenched disorder on 
the other ordered systems with a continuous degeneracy. 
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